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حکیده 


اهداف: پایش تغییرات کاربری‌ها و درک پویایی آن در یک peel Ley‏ از جایگاه 
خاصّی در مدیریت پایدار OF‏ حوضه برخوردار است. هدف تحقیق حاض استفاده از 
سنجش از دور و GIS‏ جهت igi‏ نقشة تغییرات و شناسایی انتقالات کاربری اراضی و 
پوشش زمین با به‌کارگیری ماتریس انتقال و تصاویر ماهوارة لندست در حوضة آبخیز 
دریاچۀ ارومیه می‌باشد. 

روش: جهت انجام تحقیق. از تصاویر ماهواره لندست در دور زمانی ۲۰۱۵ - ۱۹۸۸ 
استفاده گردید. بدین‌منظور پس از انجام پیش‌پردازش‌های موردنظر جهت انجام طبقه‌بندی 
از روش‌های ماشین‌بردار پشتیبان و روشیءگرا استفاده و سپس اعتبارسنجی گردیدند. 
همچنین جهت برآورد میزان انتقالات و دیگر ویژگی‌های Lege‏ آبخیز دریاچۀ ارومیه 
ابتدا ماتریس انتقالی استخراج شده و سپس طبقه‌بندی شی گرا بین دوره‌های زمانی 
۵ ارائه شد. سپس با استفاده از فرمول‌های «balay ge‏ میزان پایداری» افزایش؛ 
کاهش, تغییرات کل. تغییرات خالص و dale‏ همزمان کاربری‌های اراضی و پوشش زمین 


تاریخ دریافت: ۱۳۹۵/۰۳/۲۹ تاریخ تصویب: ۱۳۹۲/۰۶/۰۷ 
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۱۷۲ 


\. مقدمه 


مجلَة جغرافیا و توسعة ناحیه‌ای شمارة پیست و هشتم 


یافته‌ها/ نتایج: پس از ارزیابی صحت. صحت کی برای نقشه‌های حاصل از ماشین بردار 
پشتیبان و روش شی‌گرا به‌ترتیب برابر با ۹۶ و ٩۲‏ درصد و مقدار کاپای آنها به ترتیب ٩۲‏ و 
٩‏ برآورد شد که نشان‌دهندة برتری روش شی‌گرا در مقایسه با روش ماشین بردار 
پشتیبان است. در کل هر دو روش طبقه‌بندی توانستند صحت قابل‌قبولی برای نقشه‌های 
کاربری اراضی و پوشش زمین ارائه دهند. نتایج حاصل از انتقالات نشان داد به‌طور 
میانگین» ۵4 درصد از چهر؛ زمین در حوض آبخیز Seb ye‏ ارومیه در Melb‏ زمانی ۲۰۱۵-- 
۸ پایداری پوشش داشته است. که بیشترین میزان این داوم براساس مقدار این کاربری 
در فاصلة زمانی ۱۹۸۸-۲۰۱۵ مربوطه به مناطق مسکونی می‌باشد. حدود ۱٤١‏ درصد از 
سطح حوزة آبخیز دریاچۀ ارومیه به‌صورت تبادل همزمان" بوده است. همچنین سطوح آبی 
حوضة آبخیز دریاچۀ ارومیه در 8592 زمانی فوق» بیشترین ازدست‌دادگی و کمترین تبادل 
همزمان را تجربه کرده است. 

نتیجه گیری: حوضه آبخیز دریاچۀ ارومیه در این فاصلة زمانی (۱۹۸۸-۲۰۱۵) تغییرات و 
انتقالات شدیدی را تجربه کرده است» تاجایی که تنها 0A‏ درصد از چهرۀ زمین» ثابت مانده 
و قسمت‌های دیگر انواعی از انتقال‌ها را تجربه کرده‌اند. همچنین سطوح آبی و سپس 
مراتع» بیشترین آسیب‌پذیری را تجربه کرده‌اند که نشان از افزایش اراضی فاقد پوشش و 
اراضی زراعی (کشاورزی) می‌باشد. این تجزیه‌وتحلیل مارا به سنجش و تجسّم میزان 
انتقالات عمده LULC‏ درجهت برنامه‌ریزی Saal‏ حوضه آبخیز دریاچۀ ارومیه توصیه 
می کند. ۱ 

کلیدواژه‌ها: روش‌شی‌گراء ماتریس انتقال. آسیب‌پذیری و پایداری» حوضه آبخیز دریاچۀ 


ارومیه 


با توجه به نقش منابع طبیعی در زندگی بشر لازم است اطّلاع دقیقی از چگونگی وضعیت 


الاعات به‌روز, از عوامل کلیدی در برنامه‌ریزی‌هاء تصمیم‌گیری‌ها و ابزار مدیریت در هر سازمانی 


1. swap 
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وضعیت پا حالت یک شیء يا پدیده با مشاهده آن در زمان‌های مختلف است (سینگ 4 AM‏ ص. 
٩‏ در چهار dao‏ گذشته. تغییرات کاربری اراضی در ايران به یک مشکل حاد تبدیل شده است و 
اغلب تغییرات کاربری اراضی» بدون یک برنامه‌ریزی مدون و با توجه اندک به آثار زیست‌محیطی 
آنها صورت می گیرد و این باعث تشدید روند تخریب اراضی شده است. ازآنجاکه تغییرات در 
کاربری اراضی در سطوح وسیع و گسترده صورت می گیرد» تکنولوژی سنجش از راه دور یک ابزار 
ضروری و باارزش در ارزیابی تغییرات به‌دلیل پوشش مکرر و تکراری BS‏ زمین است (لوء ماسویل 
ایروندز» و Tol ge‏ ۰.۲۰۰۶ ص. OFW‏ 

به نظر می‌رسد با استفاده از تصاویر ماهواره‌ای و تکنیک‌های طبقه‌بندی آنها می‌توان اراضی‌ای را 
که دارای شباهت‌های سطحی و بازتاب مشابه هستند. در یک طبقه قرارداد و شرط asl‏ را که همان 
گروه‌بندی اراضی مشابه است. فراهم آورد (علوی‌پناه» متین‌فس و سرمدیان؛ ۱۳۸۳. صص 
۵ از نمونه تکنیک‌های طبقه‌بندی می‌توان به روش‌های ماشین بردار پشتیبان و روش 
شی گرا اشاره کرد. روش ماشین بردار پشتیبان ازجمله روش‌های غیر پارامتریک بادگیری ماشینی 
است که با تعیین یک صفح تفکیک SES‏ بهینه در فضای ویژگی داده‌های آموزشیء کلاس‌های 
مختلف را با حداکثر جدایی بین آنهاء تفکیک می‌نماید (پتروپولوس, کالایتیدز, و واردریی" ۲۰۱۰ 
صص ۱۰۷-۹۹). جذابیت ماشین‌های بردار پشتیبان با SVM‏ در قابلیت Lel‏ در حداقل‌کردن 
ریسک ساختاری یا خطای طبقه‌بندی هنگام حل dine‏ طبقه‌بندی است (تسو و مازن ۲۰۰۹. ص. 
0 در روش پیکسل پایه» واریانس و کوواریانس داده‌ها ارزیابی و فرض می‌شود که همۀ مناطق 
آموزشی از پراکنش یکسانی برخوردار هستند. درحقیقت. نمونه‌های آموزشی باید معرف کل 
کلاس‌ها باشند. بنابراین هرچه تعداد نمونه‌های آموزشی بیشتری استفاده شود تغییرهای بیشتری از 
ویژگی‌های طیفی در آن گستره قرار خواهد گرفت (علوی پناه ۳۹۰ صص ۸۳-۸۰). در مقابل» 
طبقه‌بندی شیء گرا روشی مبتنی بر قطعه‌بندی است. قطعه‌بندی تصویر, فرآیند یکپارچه‌سازی 


پیکسل‌ها براساس همکنی پدیده‌های تصویری است و براساس چهار فاکتور پهنای پنجره * تلرانس 


1. singh 

2. Lu, D., Mausel, P., Brondizio, E., & Moran, E. 
3. Petropoulos GP, Kalaitzidis C, Vadrevu KP 

4. Support vector machines 

5. Window width 
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تشابه ؛ میانگین‌وزنی "و فاکتور واریانس وزنی " کنترل می‌شود BL)‏ و استینوچر" AAAA‏ ص. AVY‏ 
قطعات باید در داخل خود همگن بوده و تنها نمایندۀ یک طبقه باشند. نه ترکیبی از چند طبقه و 
درعین‌حال باید در کل on pai‏ ناهمگنی و احتلاف بین پدیده‌های مجاور وجود داشته باشد 
(تصویربرداری Gmbh?‏ ۲۰۰۳ ص. ۲:۹). 

در دهه‌های اخیرء تغییرات سریع کاربری اراضی و پوشش زمین در حوضه درياچة ارومیه با 
پیامدهای مهمی مثل تخریب منابع طبیعی» آلودگی‌های زیستمحیطی و رشد نامناسب شهرها همراه 
بوده است. همچنین پایین‌رفتن آب دریاچة ارومیه در طی این سال‌ها می‌تواند Slee‏ از دخالت‌های 
انسانی و سیاست‌های عمومی. شرایط اقلیمی و طبیعی با به‌طورکلی انتقالات سیستماتیک و 
تصادفی کاربری اراضی و پوشش زمین باشد. آگاهی از انواع پوشش سطح زمین و نحوة استفاده از 
زمین در بخش‌های مختلف دریاچۀ ارومیه. به‌عنوان اطّلاعات پایه برای برنامه‌ریزی‌های مختلف از 
اهمّیت ویژه‌ای برخوردار است. 

تحقیق حاضر با هدف ارزیابی آماری از تغییزات کاربری اراضی و پوشش زمین در حوزه آبخیز 
دریاچۀ ارومیه طی 92 )0 زمانی ۲۷ساله انجام شد. بدین منظور از روش طبقه‌بندی ماشین بردار 
پشتیبان. شبکه‌های عصبی مصنوعی و روش شیءگرا جهت طبقه‌بندی نظارت‌شده و از ماتریس 
احتمال برای ارزیابی آماری از Bg‏ زمين در سال‌های ۱۹۸۸و ۲۰۱۵ استفاده گردیده است. 


پنطیوس, شوساس, و مکیاهیرن" (۲۰۰۶) یک تعریف آماری از تغیبرات سیستماتیک و تصادفی 
ارائه دادند. این تعریف آماری تفاوت بین انتقالات سیستماتیک و تصادفی را در تفاوت بین نسبت 
انتظاری و مشاهده‌شده از کاربری‌ها می‌داند. برایموه"(۲۰۰۳) به بررسی انتقالات سیستماتیک و 
تصادفی در شمال غنا با استفاده از تصاویر "1×١۷‏ و طبقه‌بندی نظارت‌شده پرداخت. سپس از 


ماتریس انتقال پوشش زمین با استفاده از مقايسة دو نقشه در سال‌های مختلف جهت استخراج 


. Similarity tolerance 

. Weight mean factor 

. Weight variance factor 

. Bauer, T., & Steinnocher, K. 

. Definiens Imaging Gmbh 

. Pontius, R. G., Shusas, E., & McEachern, M 
. Braimoh & et all 

. Landsat Thematic Mapper 


س یا ین ب ONDA‏ 
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نتقالات کاربری اراضی /پوشش زمین استفاده کرد. نتایج نشان داد که انتقالات تصادفی بیشتر ples‏ 
از اشغال خودبه حود زمین توسط مهاجران و گسترش گندمزارهایی جهت اسکان مجدد بوده 
همچنین تبدیل جنگل به گندمزارها تحت یک فرایند سیستماتیک از تغییرات صورت گرفته است. 
سپس آلو و پنطیوس (۲۰۰۸) به شناسایی انتقالات سیستماتیک پوشش زمین با استفاده از سنجش 
از دور و GIS‏ برای مناطق جنگلی جنوب غنا پرداختند. برای طبقه‌بندی از تصاویر لندست TM‏ و 
برای تجزیه‌وتحلیل تغییرات. از ماتریس انتقال استفاده کردند. سپس براساس انحراف‌های حاصل از 
کاربری‌ها, انتقالات سیستماتیک و راندومی را تشخیص دادند. نتایج نشان داد مناطق جنگلی 
حفاظت‌شده تحت تأثیر انتقالات سیستماتیک به زمین‌های زیر کشت رفتند. مناندهان اوده و 
پنطیوس " (۲۰۱۰) با استفاده تصاویر لندست ۱۹۸۶ و ۲۰۰۵ با به‌کارگیری ماتریس گذر به ارزیابی 
انتقالات سیستماتیک و تصادفی در منطقهٌ ولز استرالیا پرداختند. نتایج نشان داد که کمتر از ۷ درصد 
از انتقالات صورت گرفته در منطقه» به‌صورت تغییر خالص و ۲۸ درصد این تغییرات به‌صورت 
تبادل همزمان در منطقه رخ داده است. تبدیل مرتع به bly‏ بیابان به چراگاه و چراگاه به بوته‌زار از 
تبدیلات سیستماتیک در منطقه شناسایی شد. درنهایت. به‌کارگیری از ماتریس احتمال برای 
شناسایی و سنجش میزان "LUCC‏ مفید و کارآمد می‌باشد. 

پیسانت» روجر» و استامپ ۶ )۲۰۱٤(‏ در تحقیقی با استفاده از روش شیءگرا اقدام به تهية نقشة 
جنگل‌های سطح شهر کردند و نشان دادند که روش شیءگرا نقشه‌های دقیق و درعین‌حال با 
تفکیک قوی عناصر با بازتاب طیفی یکسان (مثلاً شناسایی انواع پوشش ALS‏ از هم) تولید 
می‌نماید. کرمی» خورانی» نوحه گر» شمسی. و موسوی" (۲۰۱۵) در تحقیقی با استفاده از تصاویر 
IRS‏ منطقۀ لامرد و نیز با به‌کارگیری روش شیءگرا اقدام به Ag‏ نقشه فرسايش آبکندی کردند و 
بر دقّت بالای این روش نسبت به سایر روش‌های متداول (روش ماشین بردار پشتیبان و حداکثر 
احتمال) ash‏ کردند. در Dyed‏ جدیدی نیز آرحی (۱۳۹۶) با به‌کارگیری تصاویر ماهواره‌ای 


لندست به آشکارسازی تغییرات پوشش/کاربری اراضی با استفاده از طبقه‌بندی شیمگرا پرداخت. 


1. Alo, C. A., & Pontius, R. G. 

2. Manandhar, R., Odeh, I. O., & Pontius, R. G. 

3. Land Use and Cover Change 

4. Puissant, A., Rougier, S., & Stumpf, A. 

5. Krami, A., Khorani, A., Noohegar, A., Shamsi, S. R. F., & Moosavi, V. 


Wal‏ مجلَة جغرافیا و توسعهٌ ناحیه‌ای شمارة پیست و هشتم 


نتایج نشان داد که روش طبقه‌بندی شی گرا را می‌توان در تهية نقشه پوشش زمین و تشخیص 
تغییرات استفاده کرد. 
۳ روش‌شناسی تحقیږ 
منطقة موردنظر از الگوریتم‌های مختلف استفاده گردید. این اطْلاعات ae shal‏ شامل Loss «Gain‏ 
Net Change persistence‏ و Swap‏ می‌باشند که در بخش‌های بعدی به این اصطلاحات 
پرداخته می‌شود Jade)‏ ۱ 

جدول ۱- مشخصات آماری از تغیبرات چهرة زمین 


نشانه عنوان واحد 
Gain‏ افزایش )%( 
Loss‏ کاهش )%( 
persistence‏ تداوم )%( 
Net Change‏ تغییرات خالص )%( 
Swap‏ مبادله )%( 


۳ ۲. منطقة موردمطالعه 

حوضۀ آبخیز دریاچه یک حوضة داخلی بسته است که در آن تمام رواناب های ناشی از 
بارندگی به‌سمت دریاچۀ آب‌شور مرکزی جاری می‌شوند. حوض؛ آبخیز وسیع دریاچۀ ارومیه با 
مساحت ۵۱۷۸ کیلومترمربع» مناطقی از سه استان آذربایجان‌شرقی TA)‏ درصد) آذربایجان‌غربی 
YA)‏ درصد) و کردستان (۱۰ درصد) واقع در سمت جنوب را تحت پوشش خود دارد. مختصات 
جغرافیایی این ناحیۂ وسیع ۳۰ "۳۵ تا ۳۰ ۲۸ عرض شمالی و "۰۰ 1۸ تا ۱۵ "4۶ طول شرقی 


است. شکل (۱) موقعیت این حوضه را در سطح کشور به KSE‏ استان نشان می‌دهد. 


سال پانزدهم ارزیابی آماری از انتقالات کاربری اراضی و پوشش زمین در ... \w‏ 


موقعیت حوضه آبخیز asl o‏ 


شکل ۱- موقعیت حوضة آبخیز درياچة ارومیه 
مأخذ: نگارندگان, ۱۳۹۵ 
شکل (۲) نمودار جریانی پژوهش حاضر را نشان می‌دهد. در این پژوهش تصاویر ماهوارهای 
در چندین مرحله پردازش می‌شوند. سپس با ایجاد نمونه‌های تعلیمی از منطقه در دو دوره زمانی» 
طبقه‌بندی با استفاده از روش‌های پیکسل پایه و شیگرا تهیه و ارزیابی صحت با استفاده از نقشة 
واقعیت زمینی تهیه می‌گردد. درنهایت. ماتریس انتقال تشکیل و پارامترهای آماری محاسبه خواهد 


شد. در ادامه» این مراحل به‌صورت روشن توضیح داده شده است. 


آمار و داده 
های مورد نیاز 


نقشه توپوگرافی؛ عکس 


Google «yl» 


Ea‏ بازدید زمینی و 


طبقه بندی 


استخراج آماری از ماتریس انتقال 
تغییرات کاربری اراضی و پوشش زمین 
۲۰۱۵-۸۸ 


شکل ۲- رودنمای مراحل انجام پژوهش 
مأخذ: نگارندگان» ۱۳۹۵ 


\VA 


.٤‏ مبانی نظری تحقیق 
Not‏ مواد استفاده‌شده 


lio جغرافیا و توسعة‎ álna 


شمارهة پیست و هشتم 


تصاویر سنجنده لندست از پایگاه اینترنتی 0/۳" و "USGS‏ در دوره‌های زمانی مناسب انتخاب و 


دانلود انجام شد. این تصاویر یک‌سری از تصاویر در دسترس از ماهوارۀ لندست می‌باشد که فاقد گپ 


ناشی از اصلاح کنندۂ be‏ اسکن است (حدول ۲ جهت ارزیابی آماری؛ از فرمول‌های پنطیوس 


استفاده گردید. 
گذر ردیف 
rya rya‏ 
۱۳۹۷ ۳ 
۱۳۹۷ ۳۵ 
۸ ۳ 
Yé ۸‏ 
۸ ۳۵ 
۱۹۹ ۳۳ 
۱۹۹ ۳ 
Ye 1۷‏ 
Yo ۱۳۹۷‏ 
۸ ۳۳ 
Ye ۸‏ 
yu‏ ۳۵ 
yy ۱۹۹‏ 
Ye ۱۹۹‏ 
۱۹۹ ۳۵ 


حدول ۲- مشخصات تصاویر استفاده‌شده 


۳/۳۰ 


AFAN 
7۲۱ 
VY) 
۰/۲۱ 
YYA 
2/۳۸ 


YYA 


۰10/۰ 


\0/ 
۳ 


\0/ 


۰10/ 
*۰10/ 
*10/ 


*10/ 


تاریخ تصویر 


VAM A/F 
VAME VN 
۹۸/7/۱4 
VAMA Vs 
\QAA/ Vs 


1987/06/24 
۹۸۹/۰7/۰۷ 


۲ 


۲۰ 


۲۰ 


۲۰ 


منبع دریافت تصویر 


TIF | /http://earthexplorer.usgs.gov 
TIF | /http://earthexplorer.usgs.gov 
TIF | /http://earthexplorer.usgs.gov 
TIF | /http://earthexplorer.usgs.gov 
TIF | /http://earthexplorer.usgs.gov 
TIF | /http://earthexplorer.usgs.gov 
TIF | /http://earthexplorer.usgs.gov 


tar ftp://ftp.glcf.umd.edu/glcf 
tar ftp://ftp.glcf.umd.edu/glcf 
tar ftp://ftp.glcf.umd.edu/glcf 
tar ftp://ftp.glcf.umd.edu/glcf 
tar ftp://ftp.glcf.umd.edu/glcf 
tar ftp://ftp.glcf.umd.edu/glcf 
tar ftp://ftp.glcf.umd.edu/glcf 
tar ftp://ftp.glcf.umd.edu/glcf 


1. Global land cover facility 
2. United States Geological Survey 


سال پانزدهم ارزیابی آماری از انتقالات کاربری اراضی و پوشش زمین در ... ۱۷۹ 


۶ ۲. پیش‌پردازش و موزاییک‌بندی تصویر 

برای تصحیح هندسی» تصاویر با میزان خحطای ريشة میانگین مربعات ' حدود ۰/۳۷ پیکسل ثبت 
شد و با روی‌هم‌گذاری لایه‌های Be‏ پل بر روی تصویر تطابق‌یافته» دقت بالای عمل تطابق 
هندسی تأیید شد. همچنین از روش کاهش تیرگی پدیده " که جزو خطاهای اتمسفری و خطای 
نوار نوارشدگی که از نوع خحطاهای دستگاهی می‌باشد. برای تصحیح رادیومتریک استفاده گردید. 
پدیدۀ لکه‌های ابر استفاده از تصاویر ماهواره‌ای را محدود کرده و باعث افزایش مشکل 
توف Ul oud a‏ ا EE‏ شین ما روا کش شلف ایو از 
تصاویر ماهواره‌ای در چند مرحله صورت می‌گیرد: تشخیص اب تشخیص سای اب حذف ابر و 
بارزسازی اطّلاعات (ساراینا" TE‏ صص CW‏ نتایج حاصل از این تصحیحات در شکل 
(۳) مشاهده می‌گردد. همچنین جهت موزاییک بندی از تصاویر مختصات‌دار بدون توجه به 
tas‏ ها Slee‏ گردیت. شک CE)‏ موزایکبتنی تصاویر سال sing, GUE) VA‏ 


مأخذ: نگارندگان» ۱۳۹۵ 


1. RMSE 
2. Dark Object Subtraction (DOS) 
3. Saryans 


alee ۱۸‏ جغرافیا و توسعة ناحیه‌ای شمارة پیست و هشتم 


شکل 8 موزاییک‌بندی تصاویر منطقةٌ موردمطالعه 


ag .۳ ۶‏ نمونه‌های تعلیمی. طبقه‌بندی و ارزیابی صحت 

با توجّه به بازدید میدانی و تصویر ماهوارهای گوگل‌ارث مربوط به منطقه مشاهده گردید که در 
سال ۹۸۸ شش طبقه مناطق مسکونی. جنگل» مرتع» زراعت» اراضی فاقد پوشش و سطوح آبی 
وجود داشته است. به‌منظور حذف پیکسل‌های منفرد و پراکنده در سطح تصوير طبقه‌بندی‌شده و 
همچنین به‌دست‌آوردن تصویر مطلوب و با وضوح» از تجزیه‌وتحلیل فیلتر مجوریتی" استفاده گردید. 

در این تحقیق از دو روش ماشین بردار پشتیبان" و روش شی گرا جهت طبقهبندی و از ماتریس 
انتقال» جهت شناسایی تغییرات کاربری اراضی و پوشش زمین در حوضۀ آبخیز دریاچۀ ارومیه 
استفاده گردید. همچنین برای ارزیابی صحت نقشه‌های طبقه‌بندی‌شده با قطع‌دادن (مقايسة متناظر) 
نقشه‌های طبقه‌بندی‌شده با نقشه واقعیت زمینی حاصل از مطالعات میدانی ماتریس خطا تشکیل شد 
و براساس آن» صحت کل و ضریب YS‏ محاسبه گردید. 

££ ماتریس انتقال پوشش زمین 

تجزیه‌وتحلیل آنالیز زمین با ایجاد یک ماتریس انتقال, از مقایسه دو نقشه شروع می‌شود (جدول 
AY‏ 


1. Majority Analysis 
2. Support Vector Machine(SVM) 


سال پانزدهم ارزیابی آماری از انتقالات کاربری اراضی و پوشش زمین در ... ۱۸۱ 


جدول ۳- یک ماتریس ۲ از انتقال پوشش زمین 


مأخذ: نگارندگان» ۱۳۹۵ 


در جدول بالا ردیف‌ها نشان‌دهندة 7 کلاس در سال ۱۹۸۸ است» درحالی که ستون‌ها ویژگی‌های 
کاربری اراضی در سال ۲۰۱۵ را OLS‏ می‌دهد. بنابراین علامت Cry‏ ([17) نشان‌دهندة انتقال ویژزگی 
چهرة زمین از کلاس 1 به کلاس J‏ بین سال‌های ۱۹۸۸ تا ۲۰۱۵ و Cjj‏ تداوم کلاس J‏ را در فاصلة 
سال‌های موردنظر نشان می‌دهد. همچنین علامت alee gas Cip‏ کلاس 1 در سال ۱۹۸۸ اشغال‌شده و 
جهت محاسبهٌ آن» از فرمول زیر استفاده می‌شود (پنطیوس. ۲۰۰۶): 

Ci. = Die, Ci (\) 

در این رابطه. 0 تعداد کل کلاس‌ها می‌باشد. به‌طور مشابه Ca‏ ویژگی‌ای از چهرة زمین می‌باشد 
که در سال ۲۰۱۵ ایجاد شده‌اند که به‌وسیلةٌ فرمول زیر محاسبه می‌شود: 


77 mæ 
Cj 5 سر‎ CU (Y) 


ستون Loss‏ نشان‌دهندةٌ ویژگی‌هایی از iea‏ زمین است که یک کاهش ناخالصی از کلاس 1را 
بین سال‌های ۱۹۸۸ تا ۲۰۱۵ تجربه کرده‌اند» درحالی که ردیف bagan‏ ویژگی‌هایی از چھر؟ زمین را 
نشان می‌دهد که یک رشد ناخالصی از کلاس زرا بین سال‌های ۱۹۸۸ تا ۲۰۱۵ تجربه کرده‌اند. مفهوم 
0 دلالت برافزايش و کاهش همزمان یک کلاس زمین از چشم‌انداز است Wes)‏ احیای جنگل در 
برخی مناطق و جنگل‌زدایی در مناطق دیگر). بنابراین برای محاسبة SWAP‏ نیاز به جفت‌کردن هر 
پیکسل ازدست‌رفته با پیکسل به‌دست‌آمده می‌باشد. مقدار 87720 از کلاس J‏ به‌صورت [8 می‌باشد 


که Sy pod‏ زير محاسبه می‌شود: 


۱۸۲ مجلَة جغرافیا و توسعة ناحیه‌ای شمارة پیست و هشتم 
Gy (or)‏ — وید وا — Tmin (G,‏ = 3 
اگر تغییرات خالص برابر gain) ao‏ و loss‏ مساوی باشد) باشد» Swap‏ دو برابر gain‏ و loss‏ 
۵ یافته‌های تحقیق 
۵ ۱ نتایج حاصل از طبقه‌بندی 
پس از پیش‌پردازش‌های موردنیاز و تهیۀ نمونه‌های تعلیمی, طبقه‌بندی با روش‌های ماشین بردار 
شکل )0( نشان oala‏ شده است. پس از انجام مکرر طبقه‌بندی با مقدار پارامترهای مختلف» درنهایت؛ 
بایاس» گاما تابع ES‏ پارامتر پنالتی و سطوح هرم یا پیرامید در کرنل تابع Bb‏ حلقوی به‌ترتیب برابر 
با ۵ ۰/۱2۳ ۱۰۰و « انتخاب شد. در ادامه» نقشة حاصل از کرنل حلقوی فراهم و به‌عنوان معرف 


شکل ۵ نتایج حاصل از طبقه‌بندی L‏ استفاده از روش ماشین بردار پشتیبان 
مأخذ: نگارندگان» ۱۳۹۵ 


طبقه‌بندی, انتخاب یکسری پیکسل‌های نمونه معلوم و مقایسه کلاس آنها با نتایج طبقه‌بندی لازم 
می‌باشد. بدین‌جهت تصاویر طبقه‌بندی‌شده با نمونه‌های زمینی تهیه‌شده و میزان دقت کاربری‌ها 


was)‏ کل و ضریب کاپا) به شرح جداول زیر به‌دست آمد. 


سال پانزدهم ارزیابی آماری از انتقالات کاربری اراضی و پوشش زمین در ... ۱۸۳ 


جدول -٤‏ دقت کاربری‌ها با استفاده از روش ماشین بردار پشتیبان 


پارامترهای آماری تصویر TM‏ ۱۹۸۸ تصویر ETM+‏ ۲۰۱۵ 
ضریب LW‏ ۰/۸4 ۸۰ 
دقّت کل (درصد) ae‏ ۹۱/۳۹ 


میزان تغییرات و تکرار پامترهای مؤثر برای قطعه‌بندی» جدول (۵) تهیه گردید. این جدول میزان 


می‌دهد. 
جدول ۵- فاکتورهای قطعه‌بندی 


Soom [IR 
فاکتور‎ 
نا‎ Po [oo Te va اک‎ 


پهنای پنحره 


مهم‌ترین فاکتور تأثیرگذار در قطعه‌بندی تصاویر تلرانس تشابه است. هرچه میزان این شاخص به 
۰ نزدیک‌تر گردد اندازۂ نهایی قطعات Se SF‏ می‌شود. این فرایند تا زمانی که بهترین نتیجه 
براساس نظر استفاده‌کننده به‌دست Lh‏ تکرار > dal‏ شد (شکل A‏ 


ilona ۱۸۶‏ جغرافیا و توسعة ناحیه‌ای شمارة پیست و هشتم 


T? فاد پرخ‎ 
020 40 so 120 ep 
6G-: سے‎ “i 


شکل quis‏ حاصل از طبقه‌بندی با استفاده از روش شیءگرا 


Jys‏ دقت کاربری‌ها با استفاده از روش شیءگرا 


مأخذ: نگارندگان» ۱۳۹۵ 


پارامترهای آماری تصویر TM‏ ۱۹۸۸ تصویر ETM+‏ ۲۰۱۵ 
ضریب LIS‏ ۰/۹۱ ۰/۹ 
دقت کل (درصد) AYA‏ ۹0/0 


با توجه به نتایج حاصل از جدول (۷) می‌توان گفت. روش شیءگرا به‌ترتیب با صحت کل ٩۶‏ 
و ضریب کاپای ۰/۹۲ نسبت به روش ماشین بردار پشتیبان (با صحت کل ٩۳‏ و ضریب SUS‏ 
۹ عملکرد بهتری در Le‏ نقشة کاربری اراضی و پوشش زمین داشته است: این امر به علت 
استفاده از ویژگی‌های مختلفی» ازقبیل ویژگی‌های هندسی پدیده‌ها علاوه بر ویژگی‌های طیفی آنها 


در طبقه‌بندی است. 
حدول ۷- نتایج ارزیابی صحت دو روش موردمطالعه 
مأخذ: نگارندگان» ۱۳۹۵ 
پارامترهای آماری ماشین بردار پشتیبان Seu‏ 


"AY U~ LIS ضریب‎ 


دقت کل (درصد) A.4۹۱‏ ۹۳/۹۶ 


سال پانزدهم ارزیابی آماری از انتقالات کاربری اراضی و پوشش زمین در ... ۱۸۵ 


۵ ۳ خلاصه‌سازی CVE‏ پوشش زمین 

جدول dade (A)‏ نسبت swap doss gain‏ و not change‏ برای هر پوششی از زمین 
بین سال‌های ۱۹۸۸ و ۲۰۱۵ می‌باشد. زراعت. مراتع و اراضی فاقد پوشش و شور بیشترین gain‏ 
را دارد. درحالی که مراتع» زراعت و سطوح آبی بالاترین loss‏ را دارا می‌باشند. مراتع حدود ۲۵ 
درصد. زراعت ۳ درصد و سطوح آبی ٩‏ درصد از 1088 را از چهرة زمین در منطقه تجربه کرده‌اند. 
Loss‏ در مراتع احتمالاً به‌دلیل افزایش جمعیت و افزایش زمین‌های زراعی و در سطوح آبی به‌دلیل 
حفر چاه‌ها؛ ایجاد پل و سد بر روی رودخانه‌ها و دریاچۀ ارومیه می‌باشد. بیشترین نسبت gain‏ به 
5 متعلّق به سطوح آبی است که نشان می‌دهد این کاربری» قسمت اعظمی از پوشش خود را 
ازدست داده و افزايش چندانی نداشته است. بیشترین SWAP‏ متعلّق به زراعت می‌باشد که همزمان 
با gin‏ در منطقه loss‏ نیز وجود داشته است. تغییرات در مناطق جنگلی و تاحدودی مناطق 
مسکونی تقریباً اخالص(تغییراتی که به‌صورت همزمان افزایش و کاهش کاربری اراضی را تجربه 
کرده باشد) می‌باشد. Ul‏ در کاربری‌های دیگر بیشتر به‌صورت خالص (تغییراتی که فقط افزایش یا 
کاهش پوشش زمین را تجربه کرده باشد) است. بیشترین تغیبر نسبت به مقدار (net change)‏ 
مربوط به سطوح آبی AA)‏ درصد) است. درحالی که بیشترین تغیبرات نسبت به موقعیت (SWAP)‏ 


جدول ۸ خلاصه‌سازی تغییرات چهرة زمین(درصد) 


total total > total net 
1988 2015 | 89" loss change Swap change 
مناطق مسکونی‎ 121۳۱۳ AV £/£4) 213 1۳۲ ۰۳۵۲ t/o» 
جنگل‎ ۱/۰۹ ۱/۳۷۲ ۰, ۱:۸ ۲ YA YE 
مراتع‎ 0/7.0 ۳/۸۰ VAYA | ۳ ۳۷/۳۸۱ ۳/۵ ۱ 
ely A/V ۳۷/۲۵۸ | ۲۳/۱۵۲ ۱ YV oY V/A ۱۹/۳۹۱ 
بدون پوشش‎ ۳/۳۹ ۱۳/۲۱۶ | ۱۱/۲۸۵ ۹ ۱۳/۹5۷۶ Y/OVA 4/440 
سطوح آبی‎ ۱/۷۰ VAY ۱ 44۸ 4/404 ۳ ۹/4۷ 
total ye ۱.۰ 1۰/۸۳۱ ۷۱ AVAT ۱/۶۹ 17/۱4۲ 


wI‏ مجلَة جغرافیا و توسعة ناحیه‌ای شمارة پیست و هشتم 


با توجّه به جدول (۸» هرچند که تغیبرات GE‏ در مراتع و زراعت مشهود می‌باشد. امّا به این 
نسبت تغییرات خالص این کاربری‌ها نیز زیاد می‌باشد. جنگل حالت تقریباً ثابتی را دارد. هرچندکه 
جنگل‌ها در بعضی مناطق افزایش خالص داشته» همزمان جنگل‌زدایی را هم تجربه کرده‌اند. با توجه 
به اینکه سطوح آبی افزایش خیلی کمی را در منطقه تجربه کرده است. بنابراین میزان SWAP‏ بسیار 
ضعیف و میزان تغییرات خالص و کل تقریباً برابر می‌باشد. 

.٤ ۵‏ تداوم پوشش زمین 

نسبت اختلاف پوشش‌های زمین به‌صورت ایستا بین سال‌های ۱۹۸۸ و ۲۰۱۵ در قطر جدول 
(A)‏ (به‌صورت درصد) نشان داده شده است. تقریباً ۳۱ درصد از چهرة زمین که به‌صورت مرتع 
می‌باشد» بین سال‌های ۱۹۸۸ و ۲۰۱۵ ثابت مانده است. همچنین زراعت دارای VE‏ درصد پوشش 
ثابت در این فاصلة زمانی می‌باشد. 

حدول -٩‏ ماتریس تغییرات پوشش زمین(درصد) 
deb‏ نگارندگان» ۱۳۹۵ 


۱۹۸۸/۲۰۱۵ | سطوح آبی | بدون پوشش | زراعت مراتع | جنگل | مناطق مسکونی‎ | "۲00211988 | loss 
مناطق مسکونی‎ 22 ۰/۱۰۱ | ۰/۰۵ ۶ 12 ۰ uY 1213 
جنگل‎ ۰/۹۱ ۱/۳۹۸۱ ۰ | ۹ ۰/۹ EEE VET ۸ 
مراتع‎ YYY yea | ۳۰/۸۰۲ | ۲ ۷۱ ۲ 07.0 10/0۳ 
زراعت‎ 7 ۰/۰۵ | ۰/۵۹۱۶ | ۰۹ ۱۳/۷ ۰۹ ۱۸/۷ ۳/4۰۱ 
بدون پوشش‎ er ۰/۰۰۲ | TAT ۵۹ ۱۳۹۳۹ 12۳ ۳/۳۹ ۱/۳۸۹ 
ال | ل ۱ سطوح آبی‎ | ۹ 4۹1۷ AY ۱/۷۰ ۹/۹1۸ 
total 2015 AV ۱/۳۷۲ | ۳۲/۰۸۰ ۱ | ۶ VA\Y ۱.۰ LAY) 
gain £/£4) ۰/۰۱۶ | ۱/۸۷۸ ۲ ۱۱/۸۵ ۰/۱۱ ۹/۸ 


همان‌طور که در نمودار (۱) مشاهده می‌شود مناطق جنگلی شرایط تقریباً ابتی را طی کرده‌اند 
که کمترین مقدار کاربری را در حوضهة آبخیز دریاچۀ ارومیه دارا می‌باشند. بیشترین مقدار gain‏ در 
سپس سطوح آبی می‌باشد که احتمالاً جایگزینی با زراعت و سطوح آبی داشته است. همچنین 
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نمودار ۱- نسبت افزایش. کاهش و پایداری پوشش زمین(براساس هکتار) 
مأخذ: نگارندگان» ۱۳۹۵ 


نکتۀ قابل‌ملاحظه در این نمودار وضعیت پایداری سطوح آبی می‌باشد که کاهش در منطقه 
بیشتر از پایداری این پوشش بوده است. این وضعیت نشان‌دهندة تخریب شدید در این نوع از 
چهرۂ زمین است. زراعت و مراتع شرایط تقریباً متفاوتی را تجربه کرده‌اند. به همان نسبت که 
زراعت افزایش داشته» مراتع بخش عظیمی از پوشش خود را از دست داده‌اند. 

نسبت loss‏ به persistence‏ ارزیابی آسیب‌پذیری انتقال کلاس زمین را مشخص می‌کند. 
ارزش Ip‏ (نسبت کاهش به پایداری پوشش زمین) بالاتر از یک نشان‌دهندة یک گرایش SVL‏ 
انتقال کاربری زمین به دیگر کاربری‌هایی است که مقاوم‌اند. در جدول (۱۰) Ip‏ برای کلاس سطوح 
آبی بالای یک و Lb‏ کاربری‌ها کمتر از یک می‌باشد. بنابراین این کاربری‌ها تمایل کمتری برای 
ازدست‌دادن دارند. درنتیجه» ما شاهد کاهش سطوح آبی و همچنین تاحدودی مرتع در این منطقه 
هستیم. این دو منطقه تمایل به ازدست‌دادنشان بیشتر از تمایل به افزایش در منطقه می‌باشد. 
شایانذکر است که سه کلاس مناطق مسکونی. زراعت و اراضی بدون پوشش نسبت gain-to-‏ 
(Ep) persistence‏ بالای یک دارند که نشان می‌دهد در این کاربری‌هاء Lagain‏ بیشتر به‌دلیل 
پایداری این کاربری‌ها می‌باشد. درحالی که سطوح آبی کمترین مقدار را شامل می‌شود. مقدار 
تغییرات Galle‏ به پایداری که به‌وسیلۀ فرمول Gp — ly‏ = وا[ بەد سىت le‏ برای جنگل» 
مرتع و سطوح آبی منفی می‌باشد. کاهش خالص سطوح آبی تقریباً ببعکس پایداری این چشم‌اننداز 
است. همچنین کاهش خالص زراعت با پایداری این چشم‌انداز تقریباً برابر بوده است» ولی نسبت 
به افزایش خالص نصف می‌باشد (جدول ۱۰). همان‌طور که در نمودار (Y)‏ مشاهده می‌گردد 


۱۸۸ مج جغرافیا و توسعة ناحیه‌ای شمارة بیست و هشتم 


بیشترین میزان از چهرة زمین. مربوط به تداوم این کاربری‌ها است و تبادل همزمان کاربری‌ها 


کمترین مقدار را دارا می‌باشند. 


جدول ۱۰- نسبت OP slp gp)‏ برای هر کلاس 
مأخذ: نگارندگان» ۱۳۹۵ 


gp lp np 

Yt ۱۳:۲‏ ۱/۳۷ مناطق مسکونی 
ee ۰۳۵ - 0‏ جنگل 
AD -۷‏ ۱ مراتع 
۱/۳۵ ۰/۳۷۷ ۱/۶۱ زراعت 
O/ALA M 9۱۸۰‏ بدون پوشش 
yew -VEN‏ ۲ سطوح آبی 
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swap W Absolute value of net change W persistence 


نمودار ۲- نسبت پایداری پوشش زمین با دیگر نسبت‌های جهرهة زمین (درصد) 
مأخذ: نگارندگان» ۱۳۹۵ 


با توجه به نمودار (۲)» از کل حوض آبخیز دریاچۀ ارومیه حدود ۵٩‏ درصد پایداری, ۱۶ درصد 
تبادل همزمان از تغییرات کاهشی و افزایشی و ۲۷ درصد از ger‏ زمین نیز تغییرات خالص داشته است. 
این آمارها نشان می‌دهد که در این دوره زمانی» این حوضة آبخیز تغییراتی زیادی را تجربه کرده که 
بیشتر درجهت نامناسب بوده است و همین انتقالات است که این حوضه را به یک موضوع بین‌المللی 
مبدل ساخته که با چالش‌های فراوانی روبه‌رو می‌باشد. 
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1 نتیجه گیری و پيشنهادها 

در نقشه شیءگرا» کاربری‌ها و پوشش اراضی مختلف. مانند مراتع» مناطق مسکونی» زراعت اراضی 
فاقد پوشش و سطوح آبی نسبت به دو نقشه دیگر با جزئیات بهتری طبقه‌بندی شده است؛ چراکه روش 
شیءگرا قادر است اشکال هندسی منظم. نظیر مناطق انسان‌ساخت یا اراضی زراعی را در فرایند 
سگمنت‌سازی به‌حوبی تشخیص tes‏ انتخاب پارامترها در روش‌های ماشین بردار پشتیبان اثر 
قابل‌توجهی بر صحت و زمان طبقه‌بندی دارد و ازاین‌جهت که هیچ منبع و اساسی برای انتخاب 
ogg La eal jl‏ نداشته ور bie SU‏ تجریی ple jb,‏ تونخه پیشتر می‌باشید: از ab‏ سگشت‌سازی 
روش شیءگرا نیز مستلزم انتخاب تجربی و همراه با سعی و حطای چند پارامتر است که مقدار داده‌شده 
به این پارامترها تعیین‌کنندهُ صحت و ظرافت طبقه‌بندی نهایی هستند. درمجموع. می‌توان این‌گونه 
استنباط کرد که مقايسة روش‌ها تاحدودی به منطقه و تصوير استفاده‌شده و گروه روش‌های موردمقایسه 
بستگی دارد. در اینجا هر دو روش طبقه‌بندی ماشین بردار پشتیبان و شیء گرا نشان دادند که قابلیت تهیۀ 
ec glee E‏ زا فازا مه E‏ ی شیء گرا a‏ دلیل دخالت‌دادن فاکتورهایی» نظیر 
مقیاس» شکل و بافت علاوه بر استفاده از ارزش پیکسلی, عملکرد بهتر و دقیق‌تری در طبقه‌بندی 
Bob! pie geass‏ تفا داد 

eu as‏ ما LI ere E iss‏ سا از کازاس و 
روش شیءگرا بوده است که در پژوهش حاضر نیز چنین نتیجه‌ای حاصل شد. هرچند با توجه به 
فوصت cles ya a‏ ی Leb‏ اک دو وروی clea Shall.‏ فان pled‏ 
کیفر و چیمپان ' (۲۰۰۸) قابل‌قبول و قابل‌استناد می‌باشند. امّا در پژوهش حاضر پس از مقایسة صحت 
هر دو روش. روش طبقه‌بندی شیءگرا با صحت کلی و SUIS‏ بالاتر (ضریب کاپای ۰/۹۲ دقت کل 
۷ نسبت به روش ماشین بردار پشتیبان به‌عنوان روش برتر برای منطقۀ موردهطالعه انتخاب گردید. 

در Loe‏ آبخیز دریاچۀ ارومیه با استفاده از تجزیه‌وتحلیل از ماتریس انتقال, نتایج حاصل از ارزیابی 
این ماتریس‌ها به‌ صورت زیر خحلاصه گردید. همچنین جهت نمایش این انتقالات» شکل (A)‏ 
سیگنال‌های غالب تغییرات را نشان می‌دهد. شناسابی میزان این سیکنال‌های عمده از انتقالات طبقات 


LULC‏ ما را درجهت برنامه‌ریزی آیندۀ de ge‏ آبخیز دریاچهٌ ارومیه توصیه می‌کند. 


1. Lillesand TM, Kiefer RW, Chipman JW 


alee ۱۹۰‏ جغرافیا و توسعة ناحیه‌ای شمارة پیست و هشتم 


شکل ۸ انتقالات غالب چهرة زمین ۱۹۸۸-۲۰۱۵ 


مأخذ: نگارندگان» ۱۳۹۵ 


۱ حدود ۸۱ درصد L)‏ توجه به اینکه این تغییرات جمع gain loss‏ می‌باشد که حدود ۶۱ 
درصد از سطح حوض آبخیز دریاچۀ ارومیه را دربر می‌گیرد) از چهرة زمین در حوضة آبخیز درياچة 
ارومیه تغییرات (IS‏ داشته است. از این تغیبرات» ۱۶ درصد به‌صورت SWAP‏ و مابقی تغییرات OV)‏ 
به‌صورت خالص بوده است. 

۲-در بین کاربری‌های حوضة موردنظر زراعت بیشترین SWAP‏ و مراتع بیشترین net change‏ را 
شامل بوده است. بیشترین میزان سطوح کاهشی را به‌ترتیب مراتع و سطوح آبی و بیشترین سطوح 
افزایش را به‌ترتیب زراعت و اراضی فاقد پوشش تجربه کرده‌اند. دلایل این کاهش سطوح آبی نسبت به 
مناطق دیگر می‌تواند افزایش تعداد جمعیت حفر بیش ازحد چا‌هاء ایجاد سد و حشک‌سالی‌های منطقه 
باشد. کاهش مرتع نیز می‌تواند به‌دلیل افزایش جمعیت و تبدیل کاربری مرتع به زراعت باشد. 

۳ به‌طور میانگین» ۵٩‏ درصد از چهرة زمین در حوضۀ آبخیز دریاچۀ ارومیه در فاصلة زمانی 
۵ تداوم داشته است که بیشترین مقدار از این تداوم براساس مقدار این کاربری در این فاصلة 
زمانی» مربوطه به مناطق مسکونی می‌باشد. 

٤‏ در این el‏ زمانی در حوضة آبخیز دریاچۀ ارومیه» سطوح آبی بیشترین ازدست‌دادگی و 


کمترین تبادل همزمان(5۷20) را تجربه کرده است. 
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